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A huge amount of entropy is produced at places where fresh water and seawater mix, for example 
at river mouths. This mixing process is a potentially enormous source of sustainable energy, provided 
it is harnessed properly, for instance by a cyclic charging and discharging process of porous electrodes 
immersed in salt and fresh water, respectively [D. Brogioli, Phys. Rev. Lett. 103, 058501 (2009)]. 
Here we employ a modified Poisson-Boltzmann free-energy density functional to calculate the ionic 
adsorption and desorption onto and from the charged electrodes, from which the electric work of a 
cycle is deduced. We propose optimal (most efficient) cycles for two given salt baths involving two 
canonical and two grand-canonical (dis)charging paths, in analogy to the well-known Carnot cycle 
for heat-to-work conversion from two heat baths involving two isothermal and two adiabatic paths. 
We also suggest a slightly modified cycle which can be applied in cases that the stream of fresh 
water is limited. 



Just because there's a lot of it, doesn't mean it's interesting. 
Bob Evans about water (1998). 



Charges are the invention of the devil. 
Bob Evans about ionic criticality (1998). 

I. INTRODUCTION 



Where river water meets the sea, an enormous 
amount of energy is dissipated as a result of the 
irreversible mixing of fresh and salt water. The 
dissipated energy is about 2 kJ per liter of river 
water, i.e. equivalent to a waterfall of 200m [l|. 
It is estimated that the combined power from all 
large estuaries in the world could take care of 
approximately 20% of today's worldwide energy 
demand 0]. Extracting or storing this energy 
is therefore a potentially serious option that our 
fossil-fuel burning society may have to embrace in 
order to become sustainable. However, interesting 
scientific and technical challenges are to be faced. 
So far pressure-retarded osmosis iPRO) [H3 and 
reverse electrodialysis (RED) |7Hll| have been the 
two main and best-investigated techniques in this 
field of so-called "blue energy" , or salinity-gradient 
energy. In PRO the osmotic pressure difference 
across a semi-permeable membrane is used to 
create a pressurised solution from incoming fresh 
and salt water, which is able to drive a turbine 
0-0]. In RED stacks of alternating cation- and 
anion-exchange membranes are used to generate an 
electric potential difference out of a salinity gradient 
(rl-flH. These techniques enable the generation of 
(electrical) work at the expense of the mixing of 



streams with different salinity. Actually, PRO and 
RED can be thought of as the inverse processes of 
reverse osmosis and electrodialyses, where one has 
to supply (electrical) work in order to separate an 
incoming salt-water stream in a saltier and a fresher 
stream. 



The applicability of PRO and RED are currently be- 
ing explored: a 1-2 kW prototype plant based on 
PRO was started up in 2009 in Norway [T^], and a 
5 kW RED device is planned to be upscaled to a 50 
kW demonstration project in The Netherlands [l3| . 
Interestingly, the bottleneck to large-scale applica- 
tions of both these techniques is often not the avail- 
able fuel — there is a lot of fresh and salt water — but 
rather the very large membranes that are required to 
operate at commercially interesting power outputs. 
Tailoring such membranes with a very high trans- 
port capacity and minimal efficiency losses due to 
biofouling requires advanced membrane technology. 

Recently, however, a solid-state device without mem- 
branes was constructed by Brogioli Q , who directly 
extracts energy from salinity differences using porous 
carbon electrodes immersed in an aqueous elec- 
trolyte. Due to the huge internal surface of porous 
carbon, of the order of 10^ m^ per gram of carbon, 
the capacitance of a pair of electrolyte-immersed 
porous carbon electrodes can be very large, allow- 
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ing for large amounts of ionic charge to be stored in 
the diffuse part of the double layers of the electrolytic 
medium inside the pores [l^. In fact, although the 
energy that is stored in the charged state of such 
large-area electrodes is somewhat lower than that in 
modern chargeable batteries, the power uptake and 
power delivery of these ultracapacitors is comparable 
or even larger [l^ . The capacitance of these devices 
not only scales with the contact area between the 
electrode and the electrolyte, but also with the in- 
verse distance between the electronic charge on the 
electrode and the ionic charge in the diffuse part of 
the double layer, i.e. the capacitance increases with 
the inverse of the thickness of the ionic double layer. 
As a consequence, the capacitance increases with in- 
creasing salinity, or, in other words, the potential 
increases at fixed electrode charge upon changing 
the medium from salt to fresh water. This variabil- 
ity of the capacity was used by Brogioli [3], and 
also more recently by Brogioli et al.^^, to extract 
electric work from salinity gradients without mem- 
branes. Although Sales et al. showed that the combi- 
nation of membranes and porous electrodes has some 
desirable advantages [13], we will focus here on Bro- 
gioli's experiment. 



The key concept of Ref.[lJ] is a four-stage cycle 
ABCDA of a pair of porous electrodes, together 
forming a capacitor, such that 

(AB) the two electrodes, immersed in sea water, are 
charged up from an initial state A with low 
initial charges ±Qa to a state B with higher 
charges ±Qb\ 

(BC) the salt water environment of the two elec- 
trodes is replaced by fresh water at fixed elec- 
trode charges ±Qb, thereby increasing the 
electrostatic potential of the electrodes from 

±1pB to ±1pC', 

(CD) the two highly charged electrodes, now im- 
mersed in fresh water in state C, are discharged 
back to ±Qa in state D, and finally 

(DA) the fresh water environment of the electrodes is 
replaced by salt water again, at fixed electrode 
charges ±Qa, thereby lowering the electrode 
potentials to their initial values ±ipA in state 
A. 

This cycle, during which a net transport of ions from 
salt to fresh water takes place, renders the salt water 
fresher and the fresh water saltier — although only 



infinitessimally so if the reservoir volumes are in- 
finitely large. As a consequence, the ionic entropy 
has increased after a cycle has been completed, and 
the associated free-energy reduction of the combined 
device and the two electrolyte reservoirs equals the 
electric work done by the device during the cycle, as 
we will see in more detail below. Brogioli extrapo- 
lates an energy production of 1.6 kJ per liter of fresh 
water in his device [13], equivalent to a waterfall of 
160 m, quite comparable to current membrane-based 
techniques. These figures are promising in the light 
of possible future large-scale blue-energy extraction. 
Together with the large volume of fresh and salt 
water at the river mouths of this planet, they also 
put an interesting and blessing twist to Bob Evans' 
quotes at the beginning of this article. 

Below we investigate the (free) energy and the 
performed work of electrolyte-immersed supercapac- 
itors within a simple density functional that gives 
rise to a modified Poisson-Boltzmann (PB) equation 
for the ionic double layers. By seeking analogies 
with the classic Carnot cycle for heat engines with 
their maximum efficiency to convert heat into 
mechanical work given the two temperatures of the 
heat baths, we consider modifications of Brogioli's 
cycle that may maximise the conversion efficiency 
of ionic entropy into electric work given the two 
reservoir salt concentrations. Our modification does 
not involve the trajectories AB and CD of the cycle 
where the (dis) charging electrodes are in diffusive 
contact with an electrolytic reservoir — with the 
inhomogeneously distributed salt ions "properly" 
treated grand-canonically as often advocated by 
Bob Evans [T8l - [20j . In fact, we will argue that 
the grand-canonical trajectories AB and CD at 
constant ionic chemical potential are the analogue 
of the isotherms in the Carnot cycle. Rather we 
consider to modify the constant-charge trajectories 
BC and DA (which correspond to isochores in 
a heat-engine as we will argue) by a continued 
(dis) charging process of the electrodes at a constant 
number of ions (which corresponds to an adiabatic 
(de)compression in the heat engine). In other words, 
we propose to disconnect the immersed electrodes 
from the ion reservoirs in BC and DA, treating 
the salt ions canonically while (dis)charging the 
electrodes, thereby affecting the ion adsorption and 
hence the bulk concentration from salty to fresh 
(BC) and vice versa (DA). Finally, we will consider 
a (dis) charging cycle in the (realistic) case of a finite 
volume of available fresh water, such that the ion 
exchange process renders this water brackish; the 
heat-engine analogue is a temperature rise of the 



cold bath due to the uptake of heat. 

Similar cycles were already studied theoretically 
by Biesheuvel [2l|, although not in this context of 
osmotic power but its reverse, capacitive desalina- 
tion. The "switching step" in Bicshcuvel's cycle, 
where the system switches from an electrolyte with 
a low salt concentration to an electrolyte with a 
higher salt concentration, appears to be somewhat 
different from our proposal here, e.g. without a 
direct heat-engine analogue. 
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FIG. 1: A schematic comparison of (a) the entropy- 
temperature (5, T) representation of the Carnot cycle 
of a gas that extracts mechanical work from exchanging 
heat (and thus entropy) between heat baths at two fixed 
temperatures and (b) the [N , fj,) representation (see text) 
of a cycle to extract electric work from exchanging ions 
between ionic reservoir at two fixed chemical potentials. 



II. SYSTEM AND THERMODYNAMICS 

We consider two electrodes, one carrying a charge 
Q and the other a charge —Q. The electrodes, 
which can charge and discharge by applying an ex- 
ternal electric force that transports electrons from 
one to the other, are both immersed in an aqueous 
monovalent electrolyte of volume 2V at temperature 
T. We denote the number of cations and anions 
in the volume 2V by 2N^ and 2N^, respectively. 
Global charge neutrality of the two electrodes and 
the electrolyte in the volume 2V is guaranteed if 
2iV_|- = 2N-. If the two electrodes are separated 
by a distance much larger than the Debye screening 
length — a condition that is easily met in the ex- 
periments of Ref.[l3| — then each electrode and its 
surrounding electrolyte will be separately electrically 
neutral such that Q/e = 7V_ — N+, where e is the 
proton charge and where we assume Q > without 
loss of generality. Note that this "local neutrality" 
can only be achieved provided Q/e < N+ + N^ = N, 
where the extreme case Q/e = N corresponds to an 
electrode charge that is so high that all 2N^ anions 
in the volume 2V are needed to screen the positive 
electrode and all 2N^ cations to screen the nega- 
tive one. For Q/e < N, which we assume from 
now on, we can use Q and N as independent vari- 
ables of a neutral system of the positive electrode 
immersed in an electrolyte of volume V at tempera- 
ture T, the Helmholtz free energy of which is denoted 
by F{Q, N, T, V). At fixed volume and temperature 
we can write the differential of the free energy of the 
positive electrode and its electrolyte environment as 

dF = ^diV + fdQ, (1) 

with /i — (^+-t-/i— )/2 the average of the ionic chemi- 
cal potentials fi± and ^ the electrostatic potential of 
the electrode. The last term of Eq.(IT|) is the electric 



work done on the system if the electrode charge is in- 
creased by dQ at fixed N, and hence the electrostatic 
work done by the electrode system is dW = — ^E'dQ. 
Given that is a state function, such that ^ dF — 
for any cycle, the total work done by the system dur- 
ing a (reversible) cycle equals 

W= (jidW ^- <j) *dQ = <j) fidN. (2) 

In order to be able to calculate W we thus need 
explicit cycles and the explicit equations-of-state 
^i{Q,N,T,V) and/or «'(Q, TV, T, F), for which we 
will use a simple density functional theory to be dis- 
cussed below. 

However, before performing these explicit calcula- 
tions a general statement can be made, because there 
is an interesting analogy to be made with mechan- 
ical work Wm = § pdV done by a fixed amount 
of gas at pressure p that cyclically changes its vol- 
ume and entropy (by exchanging heat). In that case 
the differential of the thermodynamic potential reads 
dU = TdS — pdV with U a state function denot- 
ing the internal energy. Since f dU = we then 
find Wm ~ § TdS. If the exchange of heat takes 
place between two heat baths at given high and low 
temperatures Tu and T/,, it is well known that the 
most-efficient cycle — the cycle that produces the 
maximum work per adsorbed amount of heat from 
the hotter bath — is the Carnot cycle with its two 
isothermal and two adiabatic (dc-)comprcssions [22j . 
If we transpose all the variables from the gas per- 
forming mechanical work to the immersed electrodes 
performing electric work, we find U <^ F ^ S ^ N , 
T ^ ^, V Q, and —p <^ where all pairs 
preserve the symmetry of being both extensive or 
both intensive. The analogue of high and low tem- 
peratures are thus high and low ionic chemical po- 
tentials fiH and fiL (corresponding to sea and river 



water, respectively), the analogue of the isothermal 
volume change is thus the (dis)charging at constant 
/X, and the analogue of an adiabatic volume change 
is (dis)charging at constant N. Therefore, the ana- 
logue of the most efficient gas cycle is the electric cy- 
cle consisting of (grand)canonical (dis)charging pro- 
cesses. Indeed, the trajectories (AB) and (CD) of 
the experimental cycle of Ref.[lJ], as discussed in 
section I, arc of a grand-canonical nature with the 
electrode in contact with a salt reservoir during the 
(dis) charging. However, the processes (BC) and 
(DA) take place at constant Q, i.e. they are equiva- 
lent to isochores in a gas cycle, instead of adiabats. 
Efficiency is thus to be gained, at least in princi- 
ple, by changing BC and DA into canonical charging 
processes. Whether this is experimentally easily im- 
plementable is, at this stage for us, an open question 
that we will not answer here. 

For the most efficient cycles, which are schemati- 
cally shown in FiglT] in the (S*, T) and the {N, /i) 
representation, we can easily calculate the work per- 
formed during a cycle. For the mechanical work of 
the gas one finds = AS" AT, with AT ^Tr-Tl 
the temperature difference and AS' the entropy that 
is exchanged between the heat baths during the 
isothermal compression and decompression. The 
analogue for the work W delivered by the electrode 
is given hy W = A/j,AiV, with A/x = hh — Ml and 
AA^ the number of exchanged ions between the reser- 
voirs during the grand-canonical (dis)charging pro- 
cesses. This result also follows directly from Eq.([5]). 
Below we will calculate AA and hence W from a 
microscopic theory. Moreover, we will also consider 
several other types of cycles. 

In the context of the thermodynamics that we dis- 
cuss here, it is also of interest to analyse the "global" 
energy flow that gives rise to the work W that the 
immersed porous electrodes deliver per (reversible) 
cycle. For this analysis it is crucial to realise that the 
device and the two salt reservoirs at chemical poten- 
tials and are considered to be at constant 
temperature T throughout, which implies that they 
are thermally coupled to a heat bath (that we call the 
"atmosphere" here for convenience) at temperature 
T. We will show that with every completed cycle, 
during which AA > ions are being transported 
from the sea to the river water, a net amount of heat 
Q > flows from the atmosphere to the two salt 
reservoirs, and that = Q in the limit that the ion 
clouds do not store potential energy due to multi- 
particle interactions. This may at first sight con- 
tradict Kelvin's statement of the Second Law ("no 
process is possible whose sole result is the complete 



£jr high salt cone'. 



low salt cone. 



FIG. 2: A scheme which illustrates the generation of 
reversible work due to a flow of AA'' particles through the 
device from a reservoir with a high salt concentration to 
a reservoir with a low salt concentration. The diflference 
in chemical potential between the two reservoirs enables 
for the extraction of an amount of reversible work W = 
(/iH — fiL)^N, which equals the decrease of the net free 
energy of the reservoirs. The increase of entropy due to 
this controlled mixing process is caused by an inflow Q = 
Qh — Ql of heat from the atmosphere, which means that 
the system tends to cool down as it performs work. Note 
that the division of this heat into two flows of different 
directions is purely suggestive, as in general there is a 
transport of heat trough the device as well. 



conversion of heat into work" [23), but one should 
realise that the cycle also involves the transport of 
ions from the sea to the river; the word "sole" in 
Kelvin's statement is thus crucial, of course. The 
analysis is based on the entropy changes ASd, ASh 
and ASl of the device, the highly-concentrated salt 
reservoir and the one with low salt concentration, 
respectively, upon the completion of a cycle. Given 
that the device returns to its initial state after a com- 
plete cycle, its entropy change vanishes and ASd = 0. 
This implies that the device, at its fixed tempera- 
ture, docs not adsorb or desorb any net amount of 
heat. During a cycle the "river" gains A A ions, and 
hence its (Hclmholtz or Gibbs) free energy changes 
by ATl = I-LlAN, while the "sea" loses AA^ ions 
such that AFh ~ —fin AN. Now the basic identity 
F = E-TS imphes that AF^ = -eAA - TASr 
and AFl = eAA - TASl, where e = E/N is the 
average energy (or enthalpy if F denotes the Gibbs 
free energy) per particle. We assume e to be in- 
dependent of density, which physically corresponds 
to the case that there are no multi-particle contri- 
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butions to the internal energy of the reservoirs, as 
is the case for hard-core systems or ions treated 
within Poisson-Boltzmann theory as ideal gases in 
a self-consistent field. The total energy in the reser- 
voirs therefore remains constant during mixing, such 
that the entropy changes of the salt reservoirs are 
TASh = ifiH - e)AN and TASl = -(ml - e)AA^- 
As a consequence of the global preservation of en- 
tropy in the reversible cycle, the ion exchange actu- 
ally drives a heat exchange whereby the sea extracts 
a net amount of heat Qh = TASh from the atmo- 
sphere, while the river dumps a net amount of heat 
Ql = —TASl into the atmosphere. Of course the 
transport of ions itself is also accompanied with a 
heat exchange in between the reservoirs, the only rel- 
evant flow is therefore the net flow of heat out of the 
atmosphere, which is Q = Q// - Ql = AfiAN = W. 
The energy flow and the particle flow of the device 
and reservoirs are tentatively illustrated in Fig. ^ 
where one should realise that the distribution of the 
heat flow from the atmosphere into the sea (Qh) and 
the river (— Ql) depends on the heat-flow from river 
to sea or vice versa, which we have not considered 
here in any detail; only the net heat flow Qh — Ql 
is fixed by global thermodynamic arguments. This 
identification of Q with W would have the interest- 
ing implication that the conversion of this work into 
heat again, e.g. by using it to power a laptop, would 
not contribute to (direct) global warming since the 
released heat has previously been taken out of the 
atmosphere [2^. It is not clear to us, however, to 
what extent this scenario is truly realistic and rele- 
vant, given that rivers, seas, and the atmosphere are 
generally not in thermal equilibrium such that other 
heat flows are to be considered. In this study we do 
not consider the heat fluxes at all, and just consider 
systems that arc small enough for the temperature 
to be fixed. 



III. MICROSCOPIC MODEL AND 
DENSITY FUNCTIONAL THEORY 



In order to calculate ii{Q,N,T,V) and 
'^{Q,N,T,V) of a charged electrode immersed 
in an electrolyte of volume V , we need a microscopic 
model of the electrode and the electrolyte. We 
consider a positively charged porous electrode with 
a total pore volume Ve, total surface area A, and 
typical pore size L. We write the total charge of the 
positive electrode as Q = e.aA with a the number 
of elementary charges per unit area. The negative 
electrode is the mirror image with an overall minus 



sign for charge and potential, see also FiglSl The 
volume of the electrolyte surrounding this electrode 
isV ~Ve + Vo, with Vo the volume of the electrolyte 
outside the electrode. The electrolyte consists of 
(i) water, viewed as a dielectric fluid with dielectric 
constant e at temperature T , (ii) an (average) 
number — [N + Q/e)/2 of anions with a charge 
— e and (iii) an (average) number iV_|_ = {N — Q/e)/2 
of cations with a charge +e. The finite pore size 
L inside the electrodes is taken into account here 
only qualitatively by regarding a geometry of two 
laterally unbounded parallel half-spaces representing 
the solid electrode, both with surface charge density 
ea, separated by a gap of thickness L filled with the 
dielectric solvent and an inhomogeneous electrolyte 
characterised by concentration profile p±{z). Here 
z is the Cartesian coordinate such that the charged 
planes are at z = and z = L. The water density 
profile Pw{z) is then, within a simple incompress- 
ibility approximation {pyj{z) -\- p+{z) -\- p-{z))v = 1 
with V a molecular volume that is equal for water 
and the ions, given hy pw[z) = 1 / v — p+{z) — P-{z). 
If the electrolyte in the gap is in diffusive contact 
with a bulk electrolyte with chemical potentials 
and /i_ of the cations and anions, we can write the 
variational grand-potential as a functional i^[p+, p-] 
given by 
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(3) 



Here the first two lines denote the ideal-gas grand 
potential of the two ionic species, with K± the ionic 
thermal wavelengths. The third line is the ideal 
water-entropy, which effectively accounts for ionic 
excluded volume interactions as it restricts the to- 
tal local ion concentration to a maximum equal 
to 1/v — of course we could have taken the much 
more accurate hard-sphere functionals closer to Bob 
Evans' heart to account for steric repulsions [25l - f29| , 
but for now we are satisfied with the more quali- 
tative lattice-gas- like description of packing j30H32 |. 
The last line of Eq.(I3|) denotes the mean- field ap- 
proximation of the electrostatic energy in terms of 
the total charge number density q{z) = P+{z) — 
P-{z) + (t{S{z) + S{z — L)) and the electrostatic po- 
tential ipiz) = kBT(f>{z)/e. Note that ip{z) is a 
functional of p±{z) through the Poisson equation 
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FIG. 3; A sketch of the two electrodes under considera- 
tion, one positively charged and the other one negatively 
charged, both in contact with an electrolyte with a com- 
pensating ionic charge. The porosity of the electrodes is 
modeled by a slit of width L filled with electrolyte in be- 
tween two solid half spaces with surface charge density 
ecr for the positive electrode and ~ea for the negative 
one, as represented by the encircled plus and minus signs. 
The ions in the gap, represented by bare plus and minus 
signs, are free to migrate throughout the pores and also 
to the volume outside the electrodes towards the other 
electrode. 



= -AnXsqiz) with As = e^/eksT the Bjer- 
rum length of water, and that '0(0) = ^ is the elec- 
trode potential. A prime denotes a derivative with 
respect to z. 

The Euler-Lagrange equations 5VL/5p±{z) = that 
describe the equilibrium concentration profiles yield 
/i-t = fc^T In (/c»sAj_/(l — 770)) : with ps the bulk 
reservoir salt concentration and 770 ~ "i-PsV the ionic 
packing fraction in the reservoir, where 4'{z) = 0. 
When the Euler-Lagrange equations are combined 
with the Poisson equation, the modified Poisson- 
Boltzmann (PB) equation with boundary conditions 
(BCs) 



4>"{z) 



k} sinh 4>{z) 



0'(z) 



1 — 770 + ^0 cosh(/)(z) ' 
— 47rAR(T; 



2 = 



0, 



(4) 

(5) 
(6) 



with = (SttAbPs)^^/^ the Debye screening 

length, is found. BC (O follows from Gauss' law on 
the surface of the electrode, and BC (jH]) from charge- 
neutrality and the symmetry with respect to the 
midplane of the gap. Note that this equation with 
accompanying BCs was already studied in Ref.jU 



'35'|. Eq.Q reduces to the standard PB equation if 
w = 0, and the large-gap case L 00 was stud- 
ied in Ref.liil. Eq. (g]) with its BCs © and (0 
forms a closed set, and once its solution is found, 
numerically in general or analytically in the special 
case that r/o = and nL — >■ 00, the required equa- 
tion of state of the electrode potential follows from 
4'(Ps,(t) = kBT(j){0)/e. Moreover, the equilibrium 
density profiles can be used to calculate the cationic 
and anionic adsorption, i.e. the excess number of 
ions per unit surface area, defined by 



z=L/2 



=0 



(P±(^) 



(7) 



Note that we integrate the profile up to z — L/2 
as required, and that our "local charge neutrality" 
assumption implies that a = T-{a,ps) — T^{a,ps)- 
Interestingly, the total surface excess of ions, defined 

by 



r(CT,ps) = r+{a,ps) + r^{(7,ps 



(8) 



is related to the total number of ions in the volume 
V^Ve + Vohy 



N = 2psV + Ar{ps,<j) 



(9) 



Below we will use expression ([7]) and ([8|) for T{ps, <t) 
to calculate N{ps,(7) from Eq.®, or to calculate 
Ps{N,a) by solving Eq.Q for ps at given N and 



Before discussing our numerical results, it is useful to 
consider the limiting case kL 3> 1 and rjQ 1^ which 
is in fact the classic Gouy-Chapman (GC) problem 
of a single, planar, charged wall in contact with a 
bulk electrolyte of point ions. In this case the PB 
equation can be solved analytically [s^ [s^l , and the 
resulting total adsorption is given by 



GC 




47r2A| 



27rAf 



2fT* 

a, tr > cr* 



(10) 



with the crossover surface charge a* = „ " . The 
crossover behavior from T oc cr^ at low cr to F = cr at 
high a signifies a qualitative change from the lin- 
ear screening regime, where the double layers ex- 
change CO- for counter ions keeping the total ion 
concentration fixed (such that T is small), to the 
nonlinear screening regime where counterion conden- 
sation takes place. For the alleged most-efficient 
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(dis)charging cycle of current interest, operating be- 
tween two ionic reservoirs with a high sah concen- 
tration ps = ph and a low one Ps ~ Pl such that 
Afj. = i^iH — l-iL = ksThi^pH / Pl)i and for which we 
argued already that the electric work per cycle reads 
W = ANAp, the GC result ^ allows for the cal- 
culation of the ionic uptake AA^ = A{r{pH,(JB) — 
r(pH, cta)) during the grand-canonical charging from 
a low charge density a a to a high one as- In 
the limit of highly charged surfaces we thus find 
AA^ = A{(7b ~ <^a), and hence the optimal work 
per unit area within the GC limit reads 

Wgc , rj., \^ Ph r „ * 

— = Kb-I- ((Jb — (Ja) tor CTA > cr . 

A PL 

(11) 

With the typical numbers ctb — a a of the order of 
nm~^, A ~ 10^ m^ per gram of porous carbon, and 
Ph/ PL = 100 one arrives at Wgc = 10 J per gram 
of carbon. Interestingly, this is substantially higher 
than Brogioli's experimental findings of only w 20 
mJ/gram per cycle. We will discuss this difference 
below. The GC limit also yields an analytic expres- 
sion for the surface potential ^gc given by 



GC 



= 2arcsinh 




With typical Debye lengths k' 
0.72 nm we find for the typical crossover surface 
charge density a* ~ 0.2 nm~^ which corresponds 
to a surface potential ~ 50 mV. For complete- 
ness we also mention the differential capacitance 
C = dQ/dvE" of an immersed electrode, which within 
the GC limit was already written by Chapman [ssj 



Cgc 
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Id* 
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47r 



cosh 



2kBT 



, (13) 



which corresponds to a two-plate capacitor with 
spacing k^^ / cosh{e^ / {2kBT)) in a dielectric 
medium characterised by its relative dielectric con- 
stant e. This result shows that the capacity indeed 
increases with the salt concentration, in agreement 
with observations that the electrode potential rises 
at fixed charge upon sweetening the surrounding 
water [ll |. 

Useful insights can be obtained from these analytic 
GC expressions. Moreover, practical linear approx- 
imations Issj . l40j and even (almost) analytical solu- 
tions |4l| - |43| for the PB equation exist in the case 



of small pore size L. Nevertheless, the pointlike na- 
ture of the ions gives rise to surface concentrations 
of counterions that easily become unphysically large, 
e.g. far beyond lOM for the parameters of interest 
here. For this reason we consider the steric effects 
through the finite ionic volume v and the finite pore 
size L below, at the expense of some numerical effort. 



IV. NUMERICAL RESULTS 



The starting point of the explicit calculations is the 
numerical solution of Eq. ^ with EC's and (O 
on a discrete z-grid of 5000 equidistant points on 
the interval z € [Q,L/2], which we checked to be 
sufficient for all values of kL that we considered. 
Throughout the remainder of this text we set v = 
with a = 0.55 nm, which restricts the total local ion 
concentration p+{z) + p~{z) to a physically reason- 
able maximum of 10 M. The Bjerrum length of water 
is set to Xb = 0.72 nm. 

We first consider a positive electrode immersed in 
a huge {V ^ Ve) ionic bath at a fixed salt concen- 
tration Ps, such that the ions can be treated grand- 
canonically. In Fig. |4] we plot (a) the electrode po- 
tential and (b) the total ion adsorption F, both as 
a function of the electrode charge number density cr, 
for three reservoir salt concentrations ps =1, 10, and 
100 mM from top to bottom, where the full curves 
represent the full theory with pore size L = 2 nm, 
the dashed curves the infinite pore limit kL ^ 1, 
and the dotted curve the analytic Gouy-Chapman 
expressions (for kL ^ 1 and v = 0) of Eas. ((TO| and 
HI]). The first observation in Fig. SJa) is that GC 
theory breaks down at surface charge densities be- 
yond le nm"^, where steric effects prevent too dense 
a packing of condensing counterions such that the ac- 
tual surface potential rises much more strongly with 
cr than the logarithmic increase of GC theory (see 
Eq. p^ ). This rise of the potential towards ~ 1 V 
may induce (unwanted) electrolysis in experiments, 
so charge densities exceeding, say, 5e nm~^ should 
perhaps be avoided. A second observation is that the 
finite pore size L hardly affects the ^'(cr) relation for 
cr > 1 nm^^, provided the steric effects are taken 
into account. The reason is that the effective screen- 
ing length is substantially smaller than L in these 
cases due to the large adsorption of counterions in 
the vicinity of the electrode. A third observation is 
that the full theory predicts, for the lower salt con- 
centrations Ps = 1 and 10 mM, a substantially larger 
5* at low cr, the more so for lower ps- This is due to 
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the finite pores size, which is not much larger than 
in these cases, such that the ionic double layers 
must be distorted: by increasing a Donnan-like po- 
tential is generated in the pore that attracts enough 
counterions to compensate for the electrode charge 
in the small available volume. Interestingly, steric 
effects do not play a large role for r{a) in Fig. IH^b), 
as the full curves of the full theory with v = are 
indistinguishable from the full theory with v — 0. 
The finite pore size appears to be more important 
for r{a), at least at first sight, at low a, where Tqc 
appears substantially lower than T from the full cal- 
culation in the finite pore. However, this is in the 
linear regime where the adsorption is so small that 
only the logarithmic scale reveals any difference; in 
the nonlinear regime at high a all curves for F coin- 
cide and hence the GC theory is accurate to describe 
the adsorption. 

We now consider the (reversible) ^'-cr cycle ABCDA 
shown in FiglSja), for an electrode with pore sizes 
L = A nm that operates between two salt reservoirs 
at high and low salt salt concentrations ps = pn = 
0.6 M (sea water) and Ps ~ Pl — 0.024 M (river wa- 
ter), respectively, such that Ap/ksT = 3.3. For sim- 
plicity we set Vo = such that the total electrolyte 
volume equals the pore volume Ve = AL/2. The 
trajectory AB represents the charging of the elec- 
trode from an initial charge density cr^ = 1 nm~^ 
to a final charge density (Jb = 2 nm~^ at ps = pn-, 
which involves an increase in the number of ions per 
unit area ATV/^ = T{aB,PH) - T{aA,PH) = 0.7 
nm^^ using Eqs.® and ([9]) which we calculate nu- 
merically with ([7]). The trajectory EC is calcu- 
lated using the fixed number of particles in state 
B, N = Nb = 2phV + AT{pH,aB), calculating a 
lower and lower value for ps for increasing cr's using 
Eq.® until ^ p^ at a ^ ac ^ 2.81 nm"^. Then 
the discharging curve CD, at fixed ps = pl is traced 
from surface charges ac down to ud = 1-97 nm~^ 
for which T{aD,PL) -V{ac,PL) = -AiV/A, i.e. the 
discharging continues until the number of expelled 
ions equals their uptake during the charging process 
AB. The final trajectory, DA, is characterised by the 
fixed number of particles in state D (which equals 
that in A), and is calculated by numerically finding 
higher and higher /Og-values from Eq.® for surface 
charges a decreasing from an to a a, where ps = Ph 
at (T = (Tyi such that the loop is closed. Note that 
all four trajectories involve numerical solutions of 
the modified Poisson-Boltzmann problem and some 
root-finding to find the state points of interest, and 
that the loop is completely characterised by pB , p^ , 
OA, and gb- Fig[5Jb) shows the concentration pro- 
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FIG. 4: (a) The surface potential ^ and (b) the ionic 
adsorption F, both as a function of the surface charge 
density a. for a porous electrode with a pore size L = 2 
nm immersed in an aqueous electrolyte of monovalent 
cations and anions at reservoir salt concentrations ps — 
{ImM, lOmM, lOOmM} from top to bottom. The la- 
bels denote results stemming from the modified Poisson- 
Boltzmann (MPB) and Gouy-Chapman theory, with the 
corresponding molecular volume v and the (finite) pore- 
size L. 



files of the anions (full curve) and cations (dashed 
curves) in the states A, B, C, and D, (i) showing 
an almost undisturbed double layer in A and B that 
reaches local charge neutrality and a reservoir con- 
centration p±{z) = Ph in the center of the pore, (ii) 
an increase of counterions at the expense of a de- 
crease of colons in going from B to C by a trade off 
with the negative electrode, accompanied by the sat- 
uration of countcrion concentration at 10 M close to 
the electrode in state C and the (almost) complete 
absence of co-ions in the low-salt states C and D, 
and (iii) the trading of counterions for colons from 
D to A at fixed overall ion concentration. 

The work done during the cycle ABCDA follows from 
either the third or the fourth term of Eq.([2]), yield- 
ing W/A = 2.3kBT nm-2 or, equivalently, W/AN = 
S.SkBT for the present set of parameters. The en- 
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FIG. 5: (a) A cycle ABCDA of the surface potential * 
and the surface charge density a of an electrode with 
pore size L = 4 nm, charging up (AB) in contact with 
an electrolyte reservoir with a high salt concentration 
Ps — pH = 0.6M from (ta = 1 nm~^ to as = 2 nm~^, 
discharging (CD) in contact with an electrolyte at low 
salt concentration ps = pL = 0.024M, while being dis- 
connected from the reservoirs (so with a fixed number of 
ions) during the additional charging (BC) and discharg- 
ing (DA). The dashed cycle ABC'D'A, for which the sur- 
face charge remains fixed upon the transfer between the 
two reservoirs, resembles the cycle of the experiments of 
Ref.[l3. (b) Counter- and co-ion concentration profiles 
P-{z) (full curves) and p+{z) (dashed curves), respec- 
tively, in the four states A, B, C, and D of (a). 



closed area of the cycle ABCDA in Fig.[S]corrcsponds 
to the amount of extracted work (up to a factor 
{eA) ), and equals the net decrease of free energy 
of the reservoirs. 

In order to compare the presently proposed type 
of cycle ABCDA with the type used in the experi- 
ments of Brogioli [iJl , where "isochores" at constant 
(7 rather than "adiabats" at constant were used to 
transit between the two salt baths, we also numer- 
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FIG. 6: The same two cycles as shown in FigOJa), but 
now in the number of ions and chemical potential (N-p) 
representation (see text), showing the constant- A'' trajec- 
tories of BC and DA in (a) and a much larger spread in 
N for the cycle ABC'D'A in (b). Whilst the solid areas 
represent the work performed during a cycle, the dashed 
areas in (b) denote losses due to the non-optimal char- 
acteristics of the cycle for the given range of A''. The 
enclosed areas give the amount of work that is extracted 
per cycle and per unit electrode area. 



ically study the dashed cycle ABC'D'A of FigElJa). 
This cycle has exactly the same trajectory AB char- 
acterised by (jO//, (TA, ctb) as before. State point C 
at PL and ac' = ctb has, however, a much smaller 
number of ions Nc = ^p^V + AT{pL,crB) than 
Nb = Nc in state B and C, because its surface 
charge ac < crc- Trajectory CD' at fixed pL is 
quite similar to CD but extends much further down 
to ctd' = <TA, where the number of ions in D' is even 
further reduced to the minimum value in the cycle 
Nqi I a = 1.0 nm~^. Finally, at fixed od' the number 
of ions increases up to A'^ by gradually increasing ps 
from PL to PH- So also this cycle is completely deter- 
mined by ph, Pl, ca, and as- The electric work W' 
done during the cycle ABC'D'A follows from Eq.(2) 
and reads W /A = 3.2kBT nm~^, which is equiva- 
lent to W'/AN' = IMbT where AA^' = A^^' - A^d' 
is the number of ions that was exchanged between 
the two reservoirs during the cycle. 

Clearly, W' > W, i.e. the Brogioli-type cycle with 
the "isocharges" BC and D'A produces more work 
than the presently proposed ABCDA cycle with 
canonical trajectories BC and DA. However, the ef- 
ficiency of ABCDA, defined as W/AN, indeed ex- 
ceeds the efficiency W'/AN' of the ABC'D'A cy- 
cle. This is also illustrated in FiglSJ where the two 
cycles ABCDA (a) and ABC'D'A (b) are shown in 
the N-p representation. Whereas the total area of 
(b) is larger than that of (a), so W' > W accord- 
ing to Eq.([2|), the larger spread in AA^' compared to 
AA renders the efficiency of (b) smaller. The work 
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3.0 
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0.55 
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TABLE I: The work W and W of cycles ABCDA and ABC'D'A, respectively, as illustrated in FigslUa) andU for 
several choices of surface charges cta and as in states A and B, for systems operating between electrolytes with high 
hand low salt concentrations pn ~ 0.6 M and pL ~ 0.024 M, for electrodes with pore size L = 4 nm. We converted 
W and W' to room temperature thermal energy units ksT, and not only express them per unit electrode area A but 
also per exchanged number of ions AA'^ and AN' during the two cycles, respectively. Note that W/ AN is a property 
of the two reservoirs, not of the charge densities of the cycle. Also note that W' / AN' depends on the volume Vo of 
electrolyte outside the electrodes, here we successively give values for the optimal situation Vo = as well as for the 
situation Vo ~ Ve. 



W is therefore less than the decrease of the free en- 
ergy of the reservoirs combined. The hatched area 
of Figinijb) denotes the work that could have been 
done with the number AA^' of exchanged ions, if a 
cycle of the type ABCDA had been used. 

The fact that W > W while W/AN > W'/AN' 
proves to be the case for all charge densities aA and 
as for which we calculated W (of an ABCDA-type 
cycle) and W (of an ABC'D'A-type cycle), at the 
same reservoirs pn and pL and the same pore size 
L as above. This is illustrated in Table 1, which 
lists W and W per unit area and per transported 
ion for several choices of a a and as- The data of 
Table 1 shows that W > by up to a factor 2, 
while W/AN > W'/AN' by up to a factor of three 
for Vo = 0, and a factor 8 for Vo = Ve- We thus 
conclude that the choice for a particular cycle to 
generate electric work depends on optimization con- 
siderations; our results show that maximum work or 
maximum efficiency do not necessarily coincide. 

Table 1 not only shows the work per area and per 
ion, but in the last column also W'/AAa with 
Act = as — cta, i.e. the work per charge that is 
put on the electrode during the charging of trajec- 
tory AB. Interestingly, in these units the work is 



V. LIMITED FRESH WATER SUPPLY 

Of course many more cycles are possible. The two 
cycles ABCDA and ABC'D'A considered so far gen- 
erate electric work out of the mixing of two very 
large reservoirs of salt and fresh water, taking up 
ions from high-salt water and releasing them in fresh 



comparable to A/.i = S.S/cbT provided aA ^ cr*, as 
also follows from Gouy-Chapman theory for highly 
charged surfaces. Note that the work per trans- 
ported charge does not equal the amount of per- 
formed work per transported ion as AN/ A is typ- 
ically much larger than Aa. Nevertheless, the fact 
that W' ~ ApA{aB — aA) gives us a handle to link 
our results with the experiments of Brogioli [l^ . 
During the experiment, the charge on the electrodes 
varies by SQ = Ae{aB — aA) ~ 0.25 mC, such that 
one arrives at an expected work of 6 pi per elec- 
trode. This agrees reasonably well with the obtained 
value of 5 pJ out of the entire system. Unfortu- 
nately, the relation between the electrostatic poten- 
tial and the charge in the experiments differs signif- 
icantly from that of our theory by at least hundreds 
of millivolts; at comparable electrostatic potentials 
the charge density in Brogioli's experiments is almost 
two orders of magnitude smaller than our theoretical 
estimates. Therefore a qualitative comparison with 
the Brogioli-cycle is at this point very hard. The rel- 
atively low experimental charge densities clarify the 
lower amount of work produced per gram of elec- 
trode, which was noted earlier in the text. Including 
the Stern layer may be a key ingredient that is miss- 
ing in the present analysis |la |. 



I 

water. Due to the large volume of the two reser- 
voirs the ionic chemical potentials pu and pL, and 
hence the bulk salt concentrations pu and in the 
reservoirs, do not change during this transfer of a 
finite number of ions during a cycle. However, there 
could be relevant cases where the power output of an 
osmo-electric device is limited by the finite inflow of 
fresh water, which then becomes brackish due to the 



11 




1 1.5 2 2.5 3 

a [nm ] 



FIG. 7; Electrode potential versus electrode charge 
density a representation of a charging-discharging cycle 
ABCA of an electrode with pore size L = 4 nm (see text) . 
The trajectory AB represents a grand-canonical charging 
process during which the system is connected to a "sea" 
with a fixed high salt concentration pn = 0.6 M, the elec- 
trolyte in the pores taking up ions. The number of ions 
during trajectory EC is kept fixed by disconnecting the 
electrodes from any ionic bath, such that further charging 
leads to more ionic adsorption at the expense of the bulk 
concentration, which has reduced to the fresh water con- 
centration pL = 0.024M in C. Trajectory CA describes 
the discharging of the electrodes in contact with a finite 
volume AV — 0.751^ of initially fresh water at concen- 
tration PL, which becomes saltier due to the uptake of 
electrode-released ions until the sea concentration of salt 
pH is reached in A. 



mixing process; usually there is enough sea water to 
ignore the opposite effect that the sea would become 
less salty because of ion drainage by a cycle. In other 
words, the volume of fresh water cannot always be 
regarded as infinitely large while the salt water reser- 
voir is still a genuine and infinitely large ion bath. 
The cycle with a limited fresh water supply is equiv- 
alent to a heat-engine that causes the temperature 
of its cold "bath" to rise due to the release of rest 
heat from a cycle, while the hot heat bath does not 
cool down due to its large volume or heat capacity. 
Here we describe and quantify a cycle ABCA that 
produces electric work by reversibly mixing a finite 
volume of fresh water with a reservoir of salt water. 

We consider a finite volume AV^ = 0.75V of fresh 
water with a low salt concentration pL = 0.024 M, 
such that the number of ions in this compartment 
equals 2plAV. This fresh water is assumed to be 
available at the beginning of a (new) cycle; its fate 
at the end of the cycle is to be as salty as the sea 
by having taken up 2{pH — pl)AV ions from the 
electrode (which received them from the sea), with 



Ph = 0.6 M the salt concentration in the sea. The 
cycle, which is represented in Fig. [71 starts with 
the electrodes connected to a large volume of sea 
water at concentration pn, charged up in state A 
at a charge density aA = 0.75 nm^^. During the 
first part AB of the cycle, the electrodes are fur- 
ther charged up until the positive one has taken up 
2((0ff — Pl)AV ions in its pores, which fixes the sur- 
face charge <tb = 2.0 nm~^ in state B. Then the 
electrodes are to be disconnected from the sea, af- 
ter which the charging proceeds in trajectory BC 
such that the increasing ion adsorption at a fixed 
total ion number reduces the salt chemical potential 
down to PL (and hence the salt concentration far 
from the electrode surface down to pl) at ac = 2.8 
nm~^ in state C. The system can then be reversibly 
coupled to the finite compartment of initially fresh 
water, after which the discharging process CA takes 
place such that the released ions cause the fresh wa- 
ter to become more salty, reaching a charge density 
a A when the salt concentration in the compartment 
of volume AV equals pn- The cycle can then be re- 
peated by replacing the compartment AV by fresh 
water again. 

The relation between the surface potential the 
charge density on the electrodes a, and the ion reser- 
voir concentration ps or the ion number N, was nu- 
merically calculated using the modified PB-equation 
(H]) with BC's dS]) and combined with the adsorp- 
tion relation ([S]), with the same parameters v = a^, 
a = 0.55 nm, L = A nm, and V = Ve = AL/2 as be- 
fore. The enclosed area in Fig. [7]gives, using Eq.(l2l), 
the net amount of (reversible) work W performed 
during a cycle, which again equals the decrease of 
the free energy of the salt water reservoir and fresh 
water volume combined. In fact, this work can be 
calculated analytically as 

W^AV [{pH - PL)kBT - pUph - ml)] , (14) 

where pn — l-iL = kBTln{pH / Pl)- This result 
agrees with the prediction by Pattle [l| for very 
small PL- For the parameters of the cycle discussed 
here, we find W/AV = 1.2 kJ per liter of fresh 
water, or W/A = 0.45fcsT nm~^. The figures show 
that the amount of work per ion that is transported 
is typically smaller than what we found for the 
Carnot-like cycle, of course. 

We may compare this reversible cycle with the 
one proposed by Biesheuvel for the reverse process, 
which is called desalination. This cycle is very sim- 
ilar to ours, except that Biesheuvel's switching step 
from sca-to river water and v. v. is actually an iso- 
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r trajectory instead of our iso~N tracjectory. This 
iso-adsorption trajectory does not seem to have a 
reversible heat engine analogue, as the degree of re- 
versibility depends on the extent to which the elec- 
trolyte can be drained out of the micropores before. 
Nevertheless, wc find agreement with the work that 
must be provided in the case of only a relatively 
small output volume of fresh water, and the expres- 
sion found by Biesheuvel exactly equals Eq. (fT4)) . 
The point we would like to stress is that irreversible 
mixing during the switching step can be prevented 
by introducing a canonical(iso-A'^) part into the cy- 
cle which enables the system to adapt to a new salt 
concentration in a time-reversible fashion, such that 
maximal efficiency is preserved. 



in the river- and sea water. However, although the 
cycles used by Brogioli, with two iso— and two 
iso— cr trajectories (where the latter are analogous 
to isochores in the heat-engine) are less efficient per 
transported ion, the total work of a "Brogioli-cycle" 
is larger, at least when comparing cycles that share 
the iso— /i charging in the sea water trajectory. We 
find, for electrode potentials ^ ~ 100 — 300 mV 
and electrode charge densities cr ~ 1 — 2 nm^'^ in 
electrolytes with salt concentrations pn = 0.6 M 
(sea water) and pL ~ 0.024 M (river water), typical 
amounts of delivered work of the order of several 
ksT per transported ion, which is equivalent to 
several fc^T per nni'^ of electrode area or several kJ 
per liter of consumed fresh water. 



VI. SUMMARY, CONCLUSION AND 
DISCUSSION 

Although substantial attempts to extract renewable 
energy from salinity gradients go back to the 1970's, 
there is considerable recent progress in this field 
stemming from the availability of high-quality mem- 
branes [11| and large-area nanoporous electrodes 
[i^ l with which economically interesting yields of 
the order of 1 kJ per liter of fresh water can be 
obtained — equivalent to a waterfall of one hundred 
meter. The key concept in the recent experiments 
of Brogioli [1J| is to cyclically (dis)charge a su- 
percapacitor composed of two porous electrodes 
immersed in sea (river) water. In this article we 
have used a relatively simple density functional, 
based on mean-field electrostatics and a lattice-gas 
type description of ionic stcric repulsions, to study 
the relation between the electrode potential ^, 
the electrode surface charge density a, the ion 
adsorption F, the ion chemical potential p, and 
the total number of ions in a (slit-like) pore of 
width L that should mimic the finite pores of the 
electrodes. With this microscopic information at 
hand, we have analysed several cycles of charging 
and discharging electrodes in sea and river water. 
By making an analogy with heat engines, for which 
the most-efficient cycle between two heat baths 
at fixed temperatures is the Carnot cycle with 
isothermal and adiabatic (de)compressions, we 
considered cycles composed of iso— /i and iso— iV 
(dis)charging processes of the electrodes. We indeed 
found that these cycles are maximally efficient in 
the sense that the work per 'consumed' ion that is 
transported from the sea to the river water during 
this cycle is optimal, given the salt concentrations 



Our calculations on the Brogioli type of cycle agree 
with experiments regarding the amount of performed 
work per cycle with respect to the variance in the 
electrode charge during (dis-) charging; each unit 
charge is responsible for an amount of work that is 
given by the difference in chemical potential between 
the two reservoirs. However, the experimental data 
concerning the electrostatic potential could not be 
mapped onto our numerical data. This could very 
well be due to the fact that the pore size in the 
experiments by Brogioli is very small such that ion 
desolvation, ion polarisability, and image charge 
effects may be determining the relation between the 
surface charge and electrostatic potential. Models 
which go beyond the present mean-field description 
are probably required for a quantitative description 
of this regime. Another ingredient in a more detailed 
description must involve the finite size of the ions 
combined with the microscopic roughness of the 
carbon. The ions in the solvent and the electrons 
(holes) in the electrode material cannot approach 
infinitely close, and the resulting charge free zone 
can be modeled by a Stern capacitance. Standard 
Gouy-Chapman-Stern (GCS) theory has successfully 
been applied to fit charge- voltage curves for porous 
carbon capacitive cells [l^, within the context of 
osmo-electrical and capacitive desalination devices. 
Extensions to GCS theory are currently being 
developed which include finite pore sizes, in order to 
obtain a physically realistic and simultaneously ac- 
curate model of the Stern layer within this geometry. 

Throughout this work we (implicitly) assumed the 
cycles to be reversible, which implies that the elec- 
trode (dis)charging is carried out sufficiently slowly 
for the ions to be in thermodynamic equilibrium 
with the instantaneous external potential imposed 
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by the electrodes. This reversibihty due to the slow- 
ness of the charging process has the advantage of 
giving rise to optimal conversion from ionic entropy 
to electric work in a given cycle. However, if one 
is interested in optimizing the power of a cycle, i.e. 
the performed work per unit time, then quasistatic 
processes arc certainly not optimal because of their 
inherent slowness. Heuristically one expects that 
the optimal power would result from the trade-off 
between reversibility (slowness) to optimize the 
work per cycle on the one hand, and fast electronic 
(dis) charging processes of the electrodes and fast 
fluid exchanges on time scales below the relaxation 
time of the ionic double layers on the other. An 
interesting issue is the diffusion of ions into (or 
out of) the porous electrode after switching on (or 
off) the electrode potential |4^|43]. Ongoing work 
in our group employs dynamic density functional 
theory |48(-[50j to find optimal-power conditions for 
the devices and cycles studied in this paper, e.g. 
focussing on the delay times between the electrode 
potential and the ionic charge cloud upon voltage 
ramps. 
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The recovery of useful energy from the otherwise def- 
inite entropy increase at estuaries, which may be rel- 
evant because our planet is so full of water, is just 
one example where one can directly build on Bob 
Evans' fundamental work on (dynamic) density func- 
tional theory, inhomogeneous liquids, electrolytes, 
interfaces, and adsorption. 
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